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Abstract 

We present near- and mid-infrared observations on the shock-cloud interaction re- 
gion in the southern part of the supernova remnant HB 21, performed with the 
InfraRed Camera (IRC) aboard AKARI satellite and the Wide InfraRed Camera 
(WIRC) at the Palomar 5 m telescope. The IRC 4 fim (N4), 7 /im (S7), and 11 /im 
(Sll) band images and the WIRC H2 v = 1 — > S(l) 2.12 fim image show similar 
diffuse features, around a shocked CO cloud. We analyzed the emission through 
comparison with the H2 line emission of several shock models. The IRC colors are 
well explained by the thermal admixture model of H2 gas — whose infinitesimal H2 
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column density has a power-law relation with the temperature T, dN ~ T dT — 
with n(H 2 ) ~ 3.9 x 10 4 cm" 3 , b ~ 4.2, and iV(H 2 ;T > 100 K) ~ 2.8 x 10 21 cm" 2 . 
We interpreted these parameters with several different pictures of the shock-cloud 
interactions — multiple planar C-shocks, bow shocks, and shocked clumps — and dis- 
cuss their weaknesses and strengths. The observed H 2 v = 1 — »■ S(l) intensity is 
four times greater than the prediction from the power-law admixture model, the 
same tendency as found in the northern part of HB 21 (Paper I). We also explored 
the limitation of the thermal admixture model with respect to the derived model 
parameters. 

Key words: HB 21, SNR 89.0+4.7, IC 443, Supernova Remnant, Infrared, Shock, 
H 2 , CO 



1 Introduction 



HB 21 (G89.0+4.7) i s a la r ge (~ 120' x 90') middle-aged (~ 5000 - 7000 



yr, lLazendic and Sland 120061 ; iByun et al.l 120061 ) supernova re mnant ( SNR) 
at a distance estimated to be from ~ 0.8 kpc to ~ 1.7 kpc ( Leahyl . 1987 



Tatematsu et al. . 1990l ; Byun et al. . 2006 ). Based on its indented, shell-like ap- 
pearance in the radio and the existence of nearby giant mol ecular clouds, it is 
thought to be interacting with a molecular cloud (cf . Fig. ED lErkes and Dickell . 



1969 



Huang and Thaddeua . Il986l : iTatematsu et al.l . Il990l ) . The first direct ev- 



idence for this interaction was the detection of broad CO emission lines near 



the edge and the center of the remnant (IKoo et all l200ll ; IByun et al.l . 120061 ) . 



The existence of such an interaction was further supported by the suggestion 
that evaporation of the cloud might be respon sible for the enhanced therma l 
X-rays seen in the central part of the remnant (ILeahy and Aschenbachl . 119961 ) . 



We performed infrared imaging observations toward two localized positions 
in HB 21, where the broad CO emission lines were observ e d (Fi g. [1]), with 
two instruments: the InfraR ed Camera (IRC , lOnaka et al.l . 120071 ) aboard a 
Japanese satellite, AKA RI (IMurakami et al.l . 120071 ) and the Wide-field In- 
fraRed Camera (WIRC, IWilson et all l2003h on the Palomar 5 m Hale tele- 
scope. From the analysis of the northern part ( "Cloud N" ) data, we found that 
the mid-infrared diffuse features originated from shocked H 2 gas, with their 
excitation conditions well described by a thermal admixture of H 2 gas, whose 
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infinitesimal H2 co lumn density has a power-law relation with the temperature 



T, dN ~ T dT ( iShinn et all 120091 . hereafter Paper I). Such H 2 excitation 
conditions are consistent with the "ankle-like" energy level population dia- 
gram (i.e. a turn-up in population for higher energies, see Fig. [7]), hitherto 
observed at the shock-cloud interaction regions (cf. § 1 of Paper I). 

Here we present the analysis of the southern portion of HB 21 ("Cloud 
S"), following the method of Paper I. The near- and mid- infrared images 
(~ 2 — 13 /mi) we obtained show diffuse features around a shocked CO cloud. 
We analyze them as emission lines of H2 gas in statistical equilibrium. We find 
the emission, as with the Cloud N case, to be well described with a power- 
law admixture model of thermal H 2 gas. We then discuss these results with 
physical pictures of the shock-cloud interaction. 



2 Observations 



We observed two specific regions (Cloud N and Cloud S in Fig. [T]), where 
slow shock s (< 20 km s~ x ) propagate into clouds of nCH.2) ~ 10 3 cm" 3 



(iKoo et all l200ll ) , using two different in struments: IRC (jOnaka et all 120071 ) 



aboard the AKARI satellite and WIRC flWilson et all 120031 ) on the Palomar 



5 m telescope. The Cloud N data were analyzed in Paper I, and the Cloud S 
data are analyzed here. Details on the observations and reduction of the IRC 
and WIRC data are described separately, below. 



2. 1 AKARI IRC observations 



AKA RI is a satellite designe d for both imaging and spectroscopy in the 



infrared ([Murakami et all 120071 ). The IRC is one of AKARPs scientific in- 
struments, which covers the wavelength range 2-30 /zm and has a ~ 10' x 10' 
field-of-view for imaging. The IRC pointed-imaging observations for Cloud 
S were performed on 2007 Jun 3rd towards (RA, Dec) = (20M6 m 07.80 s , 
+50°02'02.00") in J2000. IRC comprises three channels (NIR, MIR-S, and 
MIR-L) , each of which has three band-pass filters for imaging. Among these, 
we employed four filters from NIR and MIR-S channels for the observations; 
the MIR-L channel was not used for observing Cloud S, due to lack of observ- 
ing time. Table [1] lists the wavelength coverage and the imaging resolutions 
(r), together with pixel sizes in each channel. 

Data reduction was the same as for the Cloud N data (cf. Paper I), ex- 
cept for flat-fielding; we used a different MIR-S flat, since the dark pattern 
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seen in the channel changed around 2007 Jan 7thP~l. We obtai ned the refined 



coadd ed image through the IRC Imaging Pipeline (v. 20070104 lLorente et al. 



20071 ). Astrometric i nformation was added to the coadded images, employing 



the 2MASS catalog (ISkrutskie et al.L 120061 ) . with a matching tolerance of 1.5 



pixels. The systematic errors (~ 2 — 5%) of the calibration were included in 
the error estimation, as done for the Cloud N data. Then, for the compari- 
son between images from different bands, the pixel size was interpolated to 
1" and the spatial resolution was smoothed to ~ 7 . 43". P oint sources were 



removed applying the DAOPHOT package (IStetsonl . 119871 ) of IRAF making 
use of the simple-masking method. The final images for Cloud S are displayed 
in Figure [2j 



2.2 Palomar WIRC H2 observations 



The WIRC observations were taken together with those of Cloud N (cf. 
Paper I). We carried out the H 2 v = 1 — ► S(l) 2.12 fim narrow-band filter 
imaging observation of Cloud S, centered at (RA, Dec) = (20h:46m:23.28s, 
+49° : 54' : 16.72") in J2000, on the Palomar 5 m Hale telescope on 2005 
August 29. The WIRC is equipped with a Rockwell Science Hawaii II HgCdTe 
2K infrared focal plane array, covering a 8.7' x 8.7' field of view with a ~ 
0.25" pixel scale. The WIRC field partially covered the IRC field, due to the 
mislocation of the WIRC observation center (cf. Fig. [2]). The data reduction 
was also the same as that of Cloud N data. 50 dithered images of 30 sec 
exposure were obtained. We subtracted dark and sky background from each 
individual dithered frame and then divided by a normalized flat frame. Finally, 
the dithered frames were combined to produce the final image. 

Astrometry was obtained by matching the positions of 13 field point-sources 
with those of 2MASS catalog sources, and the positions agreed within ~ 0.3". 
Flux calibration was also done using the 2MASS catalog. We matched the 
magnitudes of 13 field point-sources with the corresponding K s magnitudes 
from the 2MASS catalog. Their correlation coefficient was 0.9956, and their 
ratio, Mwirc/Mk s , was 1.026 ± 0.011. The systematic error (~ 14%) of the 
calibration was included in the error estimation. Point sources were removed 
using DAOPHOT. The full-width-at-half-maximum (FWHM) of these sources 
was found to be ~ 1.1". 



This is described in a note at http:/ /www. ir.isas.jaxa.jp/AKARI/Observation/DataReduction/IRC/ 
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3 Results 



The final images of IRC and WI RC are displayed in Figure El together with 
a 12 C0 J = 2 -> 1 230.583 GHz (IKoo et all lioOlh image for reference. The 
peak positions of th e shocked SI and S2 clouds, where broad 12 CO J = 2 — > 1 



lines were observed (IKoo et all l200ll ). are also indicated on the images. 



3. 1 Morphology 



The IRC images (Fig. [2]) show different features from band to band, as in 
Cloud N (Paper I). The N3 and N4 images are dominated by point sources, 
while the S7 and Sll images show similar diffuse features. They do not, how- 
ever, look like bow shocks or planar shocks, unlike in Cloud N. Around the 
cloud SI, common diffuse features are seen in all the IRC images, although 
they are faint in the N3 image. However, there are no such diffuse features 
around the cloud S2. This is in contrast with Cloud N, where the shocked CO 
clouds have corresponding diffuse features in the IRC S7 and Sll bands; this 
is more interesting considering that the cloud prope rties observed fro m CO 
emission lines are similar for both clouds, SI and S2 (IKoo et all l200ll ). 



Higher extinction toward the cloud S2 than the cloud SI does not seem to 
be the reason for the absence of diffuse IRC features around the cloud S2, since 
very high column density iV(H) > 10 23 cm -2 is required for the extinction to 



be effective at ~ 10 jum (iDraind . 120031 ). Recalling that H2 emissions are the 
main source for the diffuse IRC features in the case of Cloud N (Paper I), 
this absence may be caused by the lack of H 2 gas around the cloud S2. The 
dissociation of H 2 by hot gas (> 10 6 K) can be the reason, how ever it is un 



likely since the X-ray emission is not strong around the cloud S2 (IByun et al. 



2006) . The X-ray flux (0 1 - 2 .4 keV) of HB 21 is 31.8 x 10~ 10 erg s" 1 cm 



(jLeahy and Aschenbachl . Il996l). and the clo ud S2 locates ~ 10' away from the 



central X-ray emissions (IByun et all 120061 ). which corresponds to a projected 



distance of 2 — 5 pc. With a hydrogen nuclei density of 10 4 cm -3 and an at- 
tenuating column density of 10 20 — 10 22 cm -2 , this distance corresponds to an 
effective ionization parameter, Log £ e //, ranging fr om —4 to —7, sufficie ntly 



small so that the effects of X-rays are negligible (IMaloney et all Il996l ). At 



the moment, it remains uncertain why such an absence of diffuse IRC features 
happens only to the cloud S2. 



A diffuse, looplike feature with a diameter of ~ 4' is seen in the northern 
portion of the S7 and Sll images, however, it does not seem to be related 
with the shocked CO clouds SI and S2 . The 12 CO J = 2 — > 1 map shows 
a similar looplike feature at 9.4 km s -1 (IKoo et all l200ll ). Since the looplike 
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feature looks similar in the IRC S7 and Sll bands, they may be generated 
by [Ar II] 6.99 fim and [Ne II] 12.8 /^m emission lines, which are expected to 
show similar distributions in the shock-cloud i nteraction regions considering 
their ionization potentials ( Neufeld et al. ^ 2007h; indeed, these lines haye been 



frequently observed around SNRs (e.g 



Reach et all |2002| ; iNeufeld et all 120071 ). Thermal emission from warm dust 



Arendt et al.. 1999; Oliva et aL 1999 



[~> 100 K) is another candidate. However, it seems unlikely because hot gas 
[> 10 6 K) — the heat source for the warm dust — is not abundant around the 



cloud SI (Fig. 1 of lBvun et all 12000 ) 



The WIRC H 2 v = 1 — > S(l) image only covered a small portion of 
the field of the IRC images because of the mislocation, but does include the 
cloud SI (cf. Fig. [2]). The H 2 image shows diffuse features around the cloud 
SI, similar to those seen in the IRC images. The similarity is more easily 
recognizable in Figure [3X which zooms into the area around the cloud SI. 
The RGB color image is made with N4 (blue), S7 (green), and Sll (red). 
Three elongated clumps, whose sizes are comparable with the FWHM of the 
IRC images (~ 7.4"), are apparent. Overall, their colors are red-and- yellowish, 
although the southwestern part of the features shows a little bluish color. The 
filamentary features seen in the WIRC image have a similar overall morphology 
to those seen in the RGB image. They also surround the shocked CO cloud, 
SI (cf. Fig. [2] and [3]). This geometrical relationship thus suggests that the 
diffuse infrared features seen around the cloud SI may also originate from 
shock excitation. 



3.2 Quantitative Infrared Characteristics of the Shocked Gas 



Since the cloud S2 shows no relevant feature in the IRC images and was 
not covered in the WIRC H 2 v = 1 — ► S(l) image, we analyzed the cloud SI 
only. To quantify the infrared characteristics of the cloud SI, we measured its 
intensity in the IRC and WIRC images. The regions are outlined by the two 
concentric circles (see Fig. [3D- The inner circle is the source region and the 
surrounding annular shell is the background region. To avoid any contamina- 
tion during the measu rement , possible point sources were excluded referring 



2MASS point sources (ISkrutskie et all 120061 ); they are indicated as white cir- 



cles with a red slash on Figure [3J The White circles with black shadings on the 
IRC images (Fig. [3J) are bright point sources masked out during the data re- 
duction (cf. § 12.11) . Also, the northern part of the IRC images was additionally 
excluded since the WIRC image does not fully cover this region. The masked 
area is outlined by a white tetragon with a red slash. 

Table [2] lists the measured intensities together with the IRC colors, N4/S7 
and S7/S11. The IRC intensity is the strongest in Sll, and decreases to shorter 
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wavelength. Comparing with Cloud N, the S7 and Sll intensities are greater 
by a factor of 2-3 in the cloud SI. The colors are displayed as a point in 
Figure HI The colors of the N2front in Cloud N are also displayed for compar- 
ison (cf. section 15.3.21) . The cloud SI and N2front have similar colors. The H 2 
v = 1 — > S(l) intensity was also extinction-corrected, as in Cloud N. We 
calculated the extinction factor to be ~ 0.82 (Ay = 1.8 mag), derived from 
the extinction curv e of "Milky Way, Ry = 3.1" (jWeingartner and Drainei . 
200ll ; IPraineL l2003h with the foreground hydrogen nuclei column density, 



N{ B)=N(YL l ) +gJV (H 2 )= (3.5 ± 0.4) x 10 21 cm" 2 , towards the center of HB 



21 flLee et all lipOlh 



4 Radiation Source of the Shock-Cloud Interaction Features Ob- 
served in the AKARI IRC Bands 



To interpret the infrared intensities and colors (Table [2]), we must identify 
the radiation source of the features we see in the shock-cloud interactions. 
In Paper I, based on several arguments, we concluded that shocked H 2 gas 
was the most probable explanation for the interaction features observed in 
the IRC S7, Sll, and L15 bands. In the similar manner, we here attribute the 
infrared features seen around the cloud SI in the IRC N4, S7, and Sll bands 
to shocked H 2 gas. 

Firstly, the similarity between the features seen in the IRC images and 
the Palomar WIRC H 2 d = 1 ^ S(l) image is also seen in the cloud SI 
case. Although the IRC images are rather diffuse, they definitely show three 
elongated clumpy features, which similarly locate around the cloud SI as in 
the H 2 w = 1 -> S(l) image (Fig. [3] and § !3.ip . At least, this suggests that the 
features seen in the IRC bands arise partly from H 2 emissions. Secondly, the 
same observational and theoretical arguments for H 2 emission, presented in 
Paper I, are valid over the IRC N4, S7, Sll bands (~ 3 — 13 /zm): (1) only the 
H 2 emission lines belong t o the "lines of S and H 2 ( J U p > 2)" group can produce 



spatially similar features (INeufeld et all 120071 ) ; (2) H 2 lines are the dominant 



emission from shocked molecular gas whose phy sical parameters are similar t o 
the cloud SI (v s = 20 km s~\ n(H 2 )= 10 4 cm" 3 ; iKaufman and Neufeldl[l996h . 



We also considered other possible sources for the emission, presented in Pa- 
per I: fine structure ionic lines, thermal dust continuum, Polycyclic Aromatic 
Hydrocarbons (PAHs) bands, and synchrotron radiation. Again, these do not 
seem likely either, as we discuss below. 

• Within the wavelength coverage of the IRC N4, S7, Sll bands, there are 
three strong ionic lines, [Fe II] 5.34 /im, [Ar II] 6.99 //m, and [ Ne II] 12.8 /xm , 



which have been observed in the shocked regions of SNRs (e.g. lArendt et al. 
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19991 : lOliva et all Il999l : iReach et all 120021 ; iNeufeld et all . l2007h . The ion- 



ization potentials of these ions are 7.9 eV (Fe + ), 15.8 eV ( Ar + ), and 21.6 
eV (N e + ), respectively. From their case study on four SNRs. INeufeld et al. 



( 120071 ) showed that the ions in shock-cloud interaction regions have two dis- 
tinctive spatial distributions according to their ionization potential, > 13.6 
eV and < 13.6 eV. Thus, in principle, the three ionic lines ([Fe II], [Ar II], 
and [Ne II]) can generate similar features in the IRC N4, S7, and Sll band 
images. However, the ions obs erved in the shock-c loud interaction regions 
have a low correlation with H 2 (INeufeld et all 120071 ); indeed, such low corre- 



O UCIVVCC11 |± 
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iiu. ll 2 nave 

Koo et al.. 
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1 uciitiy 

Lee et al.. 


2009) 



(e.g. 

cloud SI, which shows a good correlation between the diffuse IRC and the 
H 2 u = 1 -> S(l) features (cf. section [3TTT) . it therefore seems that ionic 
lines from [Fe II] , [Ar II] , and [Ne II] are not responsible for the features we 
measured with AKARI. 

Thermal dust continuum is not likely responsible, either. To produce a corre- 
lated features through the IRC N4, S7, and Sll bands, the dust temperature 
should be higher than 500 K. In C-shocks, though t to be opera t ing in the 
cloud SI, the dust temperature is below ~ 50 K flDraine et all Il983h . In 
addition, to heat up the dust over 500 K, there should be hot gas (> 10 6 
K), well traced in X-rays, a round the dust. Ho wever, no significant hot gas 



exists around the cloud SI (IByun et alll2006l ). 



PAHs are another candidate since they are ubiquitous and have strong , 
broad band features at 3.3, 6.2, 7.7, 8.6, 11.2, 12.7, and 16.4 /xm (jTielend . 



20081 ). However, in the similar ways presented in Paper I, they are not likely 



to be the source for the diffuse features seen in the IRC bands. PAHs are 
heated slowly and cool fast, thus it is hard t o observe the shocked PAH 
emission above the background PAH emission (jTielend . 20081): although on e 



case was claimed to detect the shocked PAH emission flTappe et all 12 006) 



hitherto, such emis sions have not been observed in shocks (Ivan Dishoecki . 



20041 : iTielensl . 120081 ) 



Synchrotron radiation is unlikely, as well, since no obvious correlation be- 
tween the 12 CO J = 2 -> 1 emission and the radio continuum was observed 
around the Cloud S (section 4.2 of lKoo et alll200ll ). 



5 Comparison to Shock Models 



To interpret the infrared intensities and colors (Tabled]), we must model 
the excitation of the features. In Paper I, we concluded that shocked H 2 gas 
was the most probable explanation, based on several arguments. Among these, 
the clearest was the similarity of the diffuse features seen in the IRC and the 
WIRC H 2 v = 1 -> S(l) images. 
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Hence, here we now analyze the band intensities as emission lines from 
shocked H 2 gas, which we calculate by applying several different shock mod- 
els. These models are the same as in Paper I, except for the non-stationary 
model. This has been excluded since some bright H 2 emission lines which fall 
into the IRC N4 band are not listed in the published models (see Table 1 of 



Flower and Pineau des Foretd . Il999l ). Further descriptions of all these models 



are given in Paper I. 



5.1 C- Shock: Isothermal H2 Gas 



It is known that the shocked H2 gas behind a planar C-t ype shock can 



be ap proximated as an isothermal and isobaric slab of gas (INeufeld et al. 



2006 f), in view of the H 2 excitation diagrams predi cted for such shocks (e.g 



Kaufman and Neufeldl . fl9~9~6l : IWikenbus et alll2000h . Hence, we first calculate 



the expected IRC colors from the emission lines of isothermal H 2 gas. Figure 
H] displays the modeled IRC colors from isothermal H 2 gas as open circles (o). 
Their trajectory moves from the lower-left corner to the upper-right corner as 
the temperature increases (i.e. becomes increasingly 'blue'). This is explained 
because pure rotational lines of H 2 , which are dominant below a few 1000 K in 
the IRC bands, have shorter wavelengths for higher upper-levels. As n(H 2 ) in- 
creases, the populations are thermalized, approaching Local Thermodynamic 
Equilibrium (LTE). 

As Figure H] shows, isothermal H 2 gas can not explain the observed IRC col- 
ors with any combination of n(H 2 ) and temperature. The ortho-to-para ratio 
(OPR) was also varied from 0.5 to 5, since the OPR is expected to be differ- 



ent from 3.0 in the inte rstellar clouds (iDalgarno et all 119731; iFlower and W att 



1984 



Lacy et al.l . ll994T ) and even in shocked gas (lTimmermannl . ll998l ; IWilgenbus et al. 



20001 ). However, these variations are not able to reproduce the observed IRC 



colors (Fig. ED- The expected IRC colors at the same temperature vary accord- 
ing to the adopted OPR; however, the locus of the IRC colors do not differ 
much from the OPR=3.0 case, that is shown in Figure [5j 



A similar result was already found for the Cloud N, and it is consistent 
with the H 2 level populations displaying an ankle-like curve (see Fig. [7]). The 
critical density of an H 2 line transition increases as the energy level of the 
upper state increases (cf. Le Bourlot et all 1999 ); hence, isothermal H 2 gas 
can only produce either a straight line (LTE) or a knee-like curve (non-LTE, 
see Fig. 1 in Paper I) in the population diagram, neither of which are observed. 
Therefore, the cloud SI also has an ankle-like H 2 population. 



This ankle-like population can be understood by the morphology of the dif- 
fuse H 2 u = 1 ^ S(l) features. These features are filamentary and surround 
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the shocked CO cloud (cf. Fig. [2] and [3]). Hence, if they are generated by shocks 
propagating into the cloud SI, a range of shock velocities are expected. Since 
the po stshock temperature of C-shock is proportional to the shock velocity as 
■u s 1,35 ( INeufeld et al.l . 120061 ) . a moderate difference in v s may result in a range 
of temperature in the shocked H2 gas. However, this explanation may not be 
valid, because the same population was also found for the N2front of Cloud N, 
whose appearance is so planar that little difference in shock speed is expected 
(Paper I). 



5.2 C -Shock: Power-law Distribution of H2 Gas Temperature 



Figure @] also displays the IRC colors calculated from the admixture model 
for H 2 gas, as filled circles (•). As can be seen, it can reproduce the observed 
IRC color ratios with an appropriate combination of n(H 2 ) and b. The derived 
parameters are listed in Table [3j iV(H 2 ; T > 100 K) is determined by scaling 
the modeled IRC intensity for the derived n(H 2 ) and b to meet the observed 
IRC intensity. The detail contributions of H 2 line emission to IRC bands for 
these parameters are listed in Table HI The "Weight" column of Table H] lists 
the weighting factor for each line to the IRC band contribution. For example, 
the Sll band intensity can be calculated as follows. 

WH 2 ) _ 0.610 7[H 2 S(2)]+ 0.921 7[H 2 ff(3)] 
MJy sr -1 10~ 4 erg s -1 cm -2 sr -1 



As Table H] shows, the pure-rotational H 2 emission lines are dominant in all 
IRC bands. Since the N4 band, in contrast to the L15 band used for Cloud N, 
is used for the color-color diagram of the SI cloud (Fig. [4J, the contribution 
from several higher-level emission lines, S(7)-S(ll), was also included, while 
that from the S(l) 17 /xm emission line was not applied (cf. Table 4 of Paper I). 
We here note that the model parameters for N2front derived from the N4/S7 
vs. S7/S11 color-color diagram (Fig. 0j), b ~ 4 and n(H 2 )~ 4 x 10 4 cm -3 , are 
both larger than those from the S7/S11 vs. S11/L15 diagram (Paper I), b ~ 3 
and n(H 2 )~ 2 x 10 3 cm~ 3 . We discuss this issue further in section IB. 3. 21 



In addition, since 12 CO y = 1 — 4.6 1 1 m em ission lines have been ob- 
served in shocked gas (e.g. iRosenthal et al.l . |2000| ) and fall into the IRC N4 
band coverage, we assessed their contribution to the band, for the power- 
law thermal admixture model with b = 4.0 and 5.0, referring the assess- 



ment of INeufeld and Yuan! (120081 ). We used the CO vibrational ene rgy state of 
Balakrishnan et al. I (l2002j ) and the CO vibrational transition rate of I Chandra et al. 



( 119961 ). We adopted the co llisional rate coefficients fo r the excita tion of CO vi- 

brational transitions by H ( Balakrishnan et al. . 20021 ) and by He ( Cecchi-Pestellini et al. 
2002 ). For excitation by H 2 , we adopted the equations [7] and [8] of iThompson 
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( 1973 ) with the parameter A of 68, the laboratory measurement of lMillikan and White 
( 119631 ). Unlike the excitation of H 2 gas, we included H as a collisional part- 
ner for the CO vibrational excitation since H excites CO vibrational levels 
(v > 0) more efficiently than He and H 2 . H excites the H 2 pure ro tational 
levels (v = 0) less efficiently than He and H 2 (ILe Bourlot et all Il999l ). hence 
including H as a collisional partner in the excitation of H 2 gas makes negligible 
effects on the predicted IRC band intensity, where the pure rotational lines 
are dominant (cf. Table H]). 



Figure M displays the results. The fractional abundance of atomic hydrogen 
to molecular hydrogen, iV(H I)/iV(H 2 ), was varied as 0, 0.01, 0.1, and 1.0, while 
that of CO was fixed as 10~ 4 . As with the H 2 vibrational states, those of CO 
also have higher collisional coefficients for a collision with atomic hydrogen 
than with He or H 2 ; hence, Figure [6] shows a sensitive dependence on the 
ratio, N(H I)/iV(H 2 ). In the range of n(H 2 )=10 3 — 10 6 cm -3 , the contribution 
to the IRC N4 band is less than 0.1, hence negligible. Furthermore, a robust 
simulation expects that 12 CO v = 1 — emission lines are much weaker than 
those of H 2 in C-shocks of preshoc k H 2 densities n(H?)= 10 4 — 10 6 cm -3 and 
shock velocities -u s =20-40 km s -1 (Kaufman and Neufeldl . Il996l ). 



The derived parameters, except the power-law index b, are a little higher 
than those previously determined towards several SNRs, where interaction 
with nearby molecular clouds is occurring. The density, n(H 2 ) = (3.9t 2 ^) xlO 4 



Koo et al.l(l200lh . The column density we derived, iV(H 2 ; T > 100 K)=(2.8lg;g) x 10 



cm , is a few times higher than the value, derived from Large Velocity Gra- 
dient (LVG) ana lysis of CO data for HB 21, of n(H 2 )= 7.0 x 10 3 cm" 3 by 

+0.2\ 
-0.5) 

cm -2 , is similarly higher than that derived towards shock-cloud interaction re- 
gions in four other SNRs (W 44, W 28, 3C 391, and IC 443), iV(H 2 )= (2.8 - 
8.9) x 10 20 cm -2 . The latter values were determined from a t wo-temperature 
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LTE fitting of pure-rotational H 2 spectra with varying OPRs (INeufeld et al. 
20071 ). Our d erived b- value of 4. 2 ln'i fa lls into the middle of the range, 3.0 



6.0, found by INeufeld and Yuan! (120081 ). These authors found the IRAC color 
ratios to be well explained with this range of power-law index b, analyzing 
Spitzer IRAC observations towards the SNR IC 443. From these three pa- 
rameters derived — n(H 2 ), iV(H 2 ;T > 100 K), and b — we also determined the 
model prediction for the H 2 v = 1 



S(l) intensity, 



;i.5ig:i)xio- 



erg s 



cm sr . It is about a factor of four smaller than the observed value (see 
Table [2] and [3]) . In contrast, for Cloud N the excess was found to be a factor 
of 17 — 33 (Paper I). We discuss these results further in §61 



Finally, we visualized the population state of the cloud SI, derived from 
the IRC color-color diagram, in Figure [7] (left). The pure-rotational levels 
which contribute to the IRC bands are designated with filled circles. Also, 
the upper level of H 2 v = 1 -> S(l) emission line is designated with a 
filled triangle; its population derived from the observed H 2 d = 1 — ■> S(l) 
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intensity, extinction corrected, is designated with a grey filled triangle with 
an error bar. The diagram shows a severe deviation of v > levels from 
v = level, it thus seems that the two temperature PTE fitting, a mode l 
usually applied for shocked H 2 gas (e.g. lRho et al.l . l200ll ; lGiannini et al.l . l2006l ). 
does not properly describe the population state of the cloud SI, even with 
varying OPR. This is caused by the low n(H.2) of the cloud SI, ~ 4 x 10 4 
' which is mu ch lower than the critic al densities for ro- vibrational H2 

19991 ). We also note that this type of 



> 10 8 cm" 3 flLe Bourlot et al. 



cm 
lines 

v > states population may not be easily recognizable with near-infrared 
ground observations, since only the population for a few lowest- J levels of 
t he v > states can be deduced due t o the atmospheric absorption (e.g. 



Burton et all Il989l ; iGiannini et all 120061 ) . Hence, for an exact derivation of 
th e level population, w e must cover full range of ~ 2 — 30 /j,m at once as 
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Rosenthal et al.l (120001 ). and space observatories are ideal and mandatory 



in this sense. 



5.3 Partially Dissociative J-shock 



As Figure [8] shows, a partially dissociating jump-shock model does not 
reproduce the observed color of the cloud SI. The observed color might appear 
to lie on a model extension to very high pressure, higher th an P = 10 1 1 cm~ 3 



K; however, this is implausible. From their CO observations. iKoo et al.l (1200 if ) 
derived n(H2) = 7.0 x 10 3 cm -3 and v s < 20km s" 1 for the cloud SI. These give 
a postshock pressure of ~ 10 s cm~ 3 K, which is more than 10 3 times lower than 
the above limit. A pressure enhancement can occur for the collision between 
mole c ular clumps and radiative shells of a remnant 



e.g. 



Moorhouse et al. 



19911 ; I Chevalier! . Il999l ); however, it is only about a factor of 20. Insufficient 



cooling time cannot solve the disagreement between the observed and modeled 
colors, either. If the postshock H 2 had not cooled as low as a few hundred K, 
then the modeled IRC colors move towards the upper-right direction in the 
color-color diagram (Fig. H]) to bluer colors. This was not observed. Overall, a 
partially dissociative J-shock does not seem to be a suitable model to explain 
the observed IRC colors. 



6 Discussion 



The observed color ratios were only reproduced by the thermal admixture 
model, as was the case for Cloud N (Paper I). Hence, we here discuss the 
derived parameters from t his model, based on pict ures for the shock-cloud 
interactions, as proposed in iNeufeld and Yuan! (l2008h and in Paper I. We also 
note here that we assumed the OPR=3.0 since no OPR information is available 
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for the cloud Si; hence, the derived parameters can be changed according to 
the adopted OPR value. 



6. 1 Nature of Molecular Shocks Seen in the Infrared 



Diffuse infrared features surround the shocked CO cloud SI (see Fig. [2] and 
[3]). For instance, the H2 n = 1 - ► S(l) image shows several filamentary 
features, as if excited by shocks propagating into the cloud core. These may 
represent distinctive planar shocks, each with different speed. From the previ- 
ous observations toward SNR molecular shock regions, the H 2 level population 
diagram has been shown to have an ankle-like curve (cf. Fig. [7]). For planar 
C-shocks to explain such populations, it generally requires two components, 
whose shock velocities are ~ 10 and ~ 3 — 50 km s" 1 , with comparable 
amounts of iV(H 2 ) (cf. iHewitt et all 120091 ) . Hence, the filamentary features 



seen in the H 2 d = 1 ^ S(l) image may originate from such a mixture 
of planar C-shocks. However, the possibility that the individual filamentary 
feature bears the ankle-like population still remains, since such a property was 
seen in a very planar filamentary feature of Cloud N (N2front). 



Neufeld and Yuan! (120081 ) showed that the h values they obtained, ~ 3.0 — 
6.0, can be explained by paraboloidal bow shocks, which are geometrical sum- 
mations of planar C-shocks (see Fig. [9]). In their picture, a paraboloidal bow 
shock, where H 2 survives the shock (i.e. T< 4,000 K), has a power-law in- 
dex b ~ 3.8. If some slower bow-shocks which do not reach 4000 K are then 
spatially averaged together, a value for b of > 3.8 is generated. 



The value b for the cloud SI was determined to be 4.2lgj, which falls into 
the range derived for bow shocks, b > 3.8. However, as discussed in Paper I, 
bow shocks should have been observed in the H 2 v = 1 —> S(l) image, if any, 
since the expected shock width for a plan ar C-shock propagati ng; into preshock 



gas of n(H 2 ) ~ 10 4 cm 3 is ~ 10 16 cm flDraine et all Il983j ). comparable to 

1.3 - 2.8) x 10 16 



the spatial resolution in the image, 

cm for the distance of ~ 0.8 — 1.7 kpc (iLeahyj . 119871 ; iTatematsu et al 



1.1" (cf. S 1272] ) 



1990 



Byun et all 120061 ). This absence of bow shock features can be explained by 
the viewing angle. If we consider the circular and filamentary appearance of 
the diffuse H 2 features around the cloud SI, it may be possible that a single 
paraboloidal bow-shock is being viewed along its symmetry axis, producing 



the circular feat ure seen in Figure 



shock picture of iNeufeld and Yuan! (120081 ). when seen face-on. 



In this case, our result fit with the bow 



This bow shock picture has a difficulty in achieving a steady state for the 
shock, however. It assumes a steady state planar C-shock at every point of 
the bow. Through the C-type shock, the n(H 2 ) can be increased up to a factor 
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of ten (e.g. Timmermann . 19981 ; Wilgenbus et al. . 2000 ). Thus, n(H 2 ) at the 
upstream of the bow would be ~ [^(H 2 ) at downstream] x 0.1 ~ [4x 10 4 ] xO.l ~ 
4 x 10 3 cm -3 (cf. Tableland Fig. [H]). Also, the preshock g as velocity into th e 



shock is known to be ~ 20 km s 1 from CO observations (IKoo et all l200ll ). 



For these preshock density and shock velocity, the time required to achieve a 
steady shoc k is known to be ~ 10 4 yr from the stu dy on the early stage of shock 
generation (IFlower and Pineau des Foretsl . Il999l ) . This time seems to be long 
for the bow shock around the cloud SI to be in a steady state, considering 
the estimated ag e of HB 21, ~ 5000 - 7000 yr flLazendic and Slanel . 120061 : 
Byun et all . 120061 ) . together with the location of the cloud SI near the edge of 
the remnant (Fig. [1]); we here note that the remnant may be older than 5000- 
7000 yr, estimat ed at the distance of 0.8 kpc, s i nce the distance is uncer tain, 
~ 0.8 - 1.7 kpc flLeahvl . Il987t iTatematsu etaD . Il990l : iBvun et all l2006h . 



In Paper I, we conjectured that a shocked clumpy interstellar medium (ISM) 
exists (cf. Fig.[H]), based on the similar b values of the N2front and N2clump re- 
gions and on the cyclodial (c uspy) feature seen in the N2clump region, together 



with numerical simulations (INakamura et all . 120061 ; IShin et all . 120081 ) . If this 



picture also holds for the SI cloud, the shocked clump must be unresolved in 
the WIRC H 2 w = 1 -> S(l) image, since the H 2 features around this cloud 
do not show any cycloidal features. However, even though this is the case, 
the absen ce of the wriggle expected for shock fronts propagating a clumpy 



!XD 

ISM (e.g. iPatnaude and Fesenl . 120051 ) still remains as an issue (cf. Fig. [3])- 



the wriggle is generated by shocks propagating further through a lower density 
medium, and vice versa. This wriggle can be unresolved in the H2 v = 1 — > 
S(l) image if its scale is small enough (< 10 16 cm). However, it is uncertain 
whether the cycloidal feature would be maintained under such a small scale. 



As noted in section I3T2"| the N4/S7 vs. S7/S11 colors of the cloud SI and 
N2front are similar (cf. Fig. H]). This is intriguing considering that they are 
physically unrelated. Their morphologies seen in H2 v = 1 —>■ S(l) images 
are also similar, i.e. filamentary, although their sizes show a few factors of 
difference (cf. Paper I and Fig. [3]). These similarities suggest that the cloud 
SI and N2front share similar shock conditions. The interstellar ultraviolet 
radiation field may contribute to this similarity; however, a more robust study 
is required. 



6.2 H 2 v = 1 — > S(l) intensity 



We estimated the H 2 d = 1 ^ S(l) intensity of the cloud SI for the 
derived model parameters — n(B.2), b, and N(B.2 ;T > 100 K) — from the mid- 
infrared IRC colors. It is about four times weaker than the observed intensity 
(Table [2] and [3]). This discrepancy is less severe than the Cloud N case, which 
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shows a factor of 17-33 difference (Paper I). However, the amount of excited 
gas, iV(H 2 ; v = 1, J = 3), required to compensate for the difference is ~ 10 14 
cm -2 in both cases (cf. Fig. [7]). 

In Paper I, we discussed two possible reasons for the discrepancy. Firstly, 
the existence of additional H 2 gas, whose temperature and density are both 
high, but whose column density is low enough to have negligible effect on the 
mid-infrared line intensities. For example, to compensate for a deficiency of 
N(v — 1, J = 3) ~ 10 14 cm -2 , we need an additional amount of H 2 gas of 
iV(H 2 ) ~ 10 16 cm -2 in LTE with T ~ 2000 K. A compact, unresolved shocked 
cloud is a candidate for such additional H 2 gas. 

The second explanation given was the omission of collisions wit h hydrogen 



atom s, which are effective in exciting the vibrational states of H 2 (cf. iNeufeld and Yuan! . 



20081 ). The cross section for excitation of H 2 by H is several orders of magnitude 
greater for rovibrational tr ansitions than it is for p ure rotational transitions 
(see Table 1 and Figure 1 in Le Bourlot et al. 1. Il999h . Hence, with only a small 



fraction of H, n(H)/n(H 2 )~ 0.025, the rovibrational transition can be dom- 
inated by collisions with H, rather than with H 2 , in the temperature range 
300-4000 K. Indeed, such a fraction of atomic gas is expe cted in interstellar 



cloud s with n(H 2 ) > 10 cm (see Table 1 and Figure 1 in lSnow and McCall 



20061 ). as well as in theoretical models for shock waves tha t are fast enough 



to pr oduce H 2 at temperatures of a few thousand K (e.g. IWilgenbus et al. 



20001 ). 



6. 3 Limitation of the Thermal Admixture Model 



6.3.1 H 2 Column Density N(H 2 ) 

In section I5T21 we mentioned that the column density iV(H 2 ; T > 100 K) 
of the cloud SI is a few times higher than those of other SNRs (W 44, W 28, 
3C 391, and IC 443). The former is iV(H 2 ; T > 10 0K) = (2.8l^) xl0 21 cm" 2 , 



while the latter are iV(H 2 ) = (2.8-8.9) x 10 20 cm" 2 (INeufeld et all 120071 ). This 
seems unreasonable considering that the n(H 2 ) of t he cloud SI, (3.9 l 2 9)xl0 4 



cm" 3 , is lower than that of IC 443, ~ 10 7 cm" 3 (INeufeld and Yuanl . 120081 ). 
However, there is one important point to note before the comparison: those 
H 2 column densities are expectations estimated using different methods. 

The total column density of H 2 gas at a few 100 K is mainly determined 
by the values ofi; = 0J = 0,l levels, since other levels have much lower 
populations (cf. Fig. [7J. However, we cannot obtain the column densities of 
these levels directly from emission lines, because no transition to a lower state 
is allowed; therefore, we must estimate the column densities of J = 0, 1 levels 
from the observed column densities of J > 1 levels. The different methods 
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used for this estimation causes the discrepancy between the cloud SI and IC 
443, mentioned above, as we discuss below. 



The iV(H 2 ; T > 100 K) of the cloud SI, (2.8lftg) x 10 21 cm" 2 , is estimated by 
applying the thermal admixture model to the observed IRC intensities, while 
the iV(H 2 ) of IC 443, ~ 5.0 x 10 20 c m" 2 , was estimated w ith two temperature 
LTE fitting with varying the OPR (INeufeld et all 120071 ). If we estimate the 
iV(n 2 ) of the cloud SI using the same two-temperature LTE fitting, we would 
obtain a lower values. Figure [7] (top-right) displays the result of such fitting. 
The fitting was applied only to the H 2 levels whose emission lines contribute 
to the IRC bands (cf. Table HJ), and returned a column density of ~ 7.3 x 10 19 
cm" 2 . This is about seven times less than IC 443, and now does not cause the 
discrepancy mentioned above. 



Also, this column density is about 4-0 times smaller than the estimation 
using the thermal admixture model. This difference stems from the fact that 
the thermal admixture model does not estimate the column densities of J = 
0, 1 with a linear extrapolation in the population diagram as the LTE fitting 
does (cf. Fig. [7]); it estimates the column densities with a curved population, 
defined by the equation dN ~ T~ b dT (100 K < T < 4000 K). Such a difference 
is smaller when longer wavelength IRC bands are used for the column density 
estimation. Figure [7J (bottom-right) shows the results of the two temperature 
LTE fitting for the population of the N2front, determined using the IRC S7, 
Sll, and L15 bands (Paper I); the result shows a column density of ~ 7.5 x 10 19 
cm -2 , which is about 3.3 times smaller than the estimation using the thermal 
admixture model, ~ 2.5 x 10 20 cm -2 (Paper I). This trend is reasonable since 
the longer wavelength IRC bands determine the level populations of lower-J 
states in v = 0, which results in a less severe extrapolation for the column 
densities, N(J = 0, 1). 



6.3.2 H 2 Density n(H 2 ) and power-law index b 

To compare the N4/S7 and S7/S11 colors of the cloud SI with those of 
N2front in Cloud N (Fig.HJ), we further determined the N4 intensity of N2front, 
which was not measured in Paper I because of strong point source contami- 
nation. In order to remove this contamination, we additionally masked point- 
source dominated areas seen in the N4 band. Table [5] shows the measured 
intensities. Since some areas are additionally masked unlike Paper I, the in- 
tensities of S7, Sll, and L15 are a little different from those listed in Paper I; 
however, the colors S7/S11 and S11/L15 are almost unchanged (cf. Fig. 7 in 
Paper I and Fig. [TU1) . Therefore, we think the point-source masking was done 
properly. 

Figure [TU] shows that the model parameters, b and n(H 2 ), obtained for 
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N2front depends on which IRC bands are used for the color-color diagram. The 
diagram of shorter wavelength bands (N4, S7, Sll; Fig. [TOl leftmost) returns 
larger-b and larger-n(H 2 ) values than the diagram of longer wavelength bands 
(S7, Sll, L15; Fig. M rightmost), and the diagram of N4/S7 vs. S11/L15 
(Fig. [10] middle) returns the middle values of the former two diagrams'. This 
inconsistency may be caused by the intrinsic property of shocked H2 gas. In 
other words, the whole level population of shocked H2 gas may not be fully 
described by the power-law thermal admixture model with only one set of b 
and n(H 2 ). 

To check this possibility, in Figure [JUJ we overplotted the IRC colors of 
Orion Molecular Cloud-1 (OMC-1), whe re extensive emis s ion li nes of shocked 



H2 gas were observed over 2.5 — 30 /iin ( [Rosenthal et all 120001 ) . The OMC-1 
emissions were adjusted to experience the same extinction with HB 21, to be 
placed in the model grid for N2front. Interestingly, OMC-1 also shows the 
same trend for the model parameters, b and n(H 2 ), as the case for N2front. 
This suggests that the variations of b and n(H2) are needed for the thermal 
admixture model to describe the whole level population of shocked H2 gas. 
No significant variation of b was s een in the SNR IC 443, where the thermal 



admixture model was applied first (INeufeld and Yuanl . l2008l ). It may be caused 



by the narrow wavelength coverage of the bands they used (Spitzer IRAC; 
~ 3 — 8 /im), which missed the longer wavelength information we used. We 
here note that the model parameters — b and n(H2) — must be obtained from 
the same band images for their comparison between different shocked regions, 
since the parameters are likely dependent on the wavelength. 



7 Conclusion 



We have observed a shock-cloud interaction region in the SNR HB 21 at 
near- and mid-infrared wavelengths, with the WIRC at the Palomar telescope 
and the IRC aboard the AKARI satellite. The IRC N4, S7, and Sll band 
images and the WIRC H 2 v = 1 — > S(l) image reveal similar diffuse features, 
which surround the shocked CO cloud SI. However, there are no infrared 
diffuse features seen around another shocked CO cloud S2. Lack of shocked 
H 2 gas may cause this absence, but why it happens only for the cloud S2 is 
uncertain. 

We found that the IRC colors of the cloud SI are well explained by an 
admixture model of H 2 gas temperatures, whose infinitesimal column density 
varies as dN ~ T~~ h dT. Three physical parameters — n(H.2), b, and iV(H2 ; T > 
100 K) — were derived from this thermal admixture model (cf. Table [3]). These 
can be understood with multiple planar C-shocks whose velocities are different. 
Alternatively, the derived b value (~ 4.2) can be understood through a bow 
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shock picture, if we are looking at a single bow shock along the symmetry axis. 
However, this picture has a difficulty in achieving a steady state. A shocked 
clumpy ISM picture, conjectured in Paper I, remains as a possible explanation, 
but the absence of the wriggle, expected for shock fronts propagating a clumpy 
medium, in the filamentary features seen in H2 v — 1 — ► S(l) image remains 
as an issue. The model parameters, b and ri(H 2 ), obtained for the cloud SI 
and N2front (cf. Fig. H]) are very similar, which means that these clouds share 
similar shock conditions. 

We also compared the observed H 2 w = 1 -> S(l) intensity to that 
predicted from the power-law admixture model. It is about four times greater. 
This excess might be caused by either an additional component of hot, dense 
H2 gas (which has low total column density), or by the omission of collisions 
with hydrogen atoms in the power-law admixture model (which results in an 
under-prediction of the near-IR line intensity). 

The limitation of the thermal admixture model is explored with respect 
to the derived model parameters. The iV(H 2 ) estimation of the model shows 
a smaller difference with those of two temperature LTE fitting, when longer 
wavelength IRC bands are used for the determination of model parameters. 
Investigating the infrared colors of N2front and OMC-1 in the four IRC bands 
(N4, S7, Sll, and L15), we found that the thermal admixture model cannot 
describe the whole H 2 level population with only one set of b and n(H 2 ); 
the shorter wavelength bands returns higher-6 and higher-n(H 2 ). This tells 
we must use the same bands in determining the model parameters, for the 
comparisons of the shocked H 2 gas' properties. 
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Fig. 1. Regions observed by AKARI. The IRC pointing-observation regions ("Cloud 
N" and "Cloud S") are overlaid as two boxes (~ 10' x 10') on the 1420 MHz 
radio continuum image of HB 21, obtained by using the synthesis telescope at the 
Dominian Radio Astrophysical Observatory. The Palomar WIRC observations were 
performed toward a similar reg ion. In this paper, we present the Cloud S data. The 
Cloud N data were presented in IShinn et all (|2009h . The 1420 MHz radio continuum 
image is kindly provided by T. L. Landecker. 
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Fig. 2. The AKARI IRC and Palomar WIRC images of Cloud S. See Table CD for 
the band definitions for the IRC images, {upper-panels) IRC N3, N4, a nd S7 band 
imag es, {lower-panels) IRC Sll band, 12 C0 J = 2 -» 1 230.583 GHz dKoo et all 
200 ll ). and H2 v = 1 — > S(l) 2.122 /im images. The peak positions ("SI" and 



"S2"), where b road CO molecul ar lines were observed, are indicated with a '+' over 
all images (cf. iKoo et al.l . l200ll ). Bright point-sources were masked out, shown by 
the white circles in the IRC images. The center position of the WIRC H2 v = 1 — > 
S(l) image was mislocated, hence the partial image is shown here. A scale bar at 
the upper- left corner shows 1'. 
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Fig. 3. The IRC RGB image (left) and the WIRC H 2 v = 1 -> S(l) image (right) 
of the cloud SI. The RGB image is composed of Sll (R), S7 (G), and N4 (B) band 
images, i.e. 11 ^m + 7 /xm + 4 /iin. All colors scale linearly, and fully cover the 
dynamic range of the diffuse features. The annular regions selected for the intensity 
measurement are indicated as two concentric white circles. The inner circle is the 
source region and the outer annular region is the background region. Circular areas 
around possible point sources were excluded during the intensity measurement to 
avoid possible contamination. These areas are indicated as white circles with a red 
slash. Bright point-sources are masked out, and their positions are indicated by 
white circles with black shading. For the comparison between the IRC and WIRC 
intensity, a tetragonal region is excluded from the IRC images. The cross ('+') 
indicates the peak positions ("SI"), where broad CO molecular lines were observed, 
as in Figure El 
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Fig. 4. The IRC color-color diagram for the cloud SI, with the colors of the N2front 
in Cloud N plotted for comparison (cf. section I6.3.2j) . The axes represent the ratio 
of the intensities in the corresponding IRC bands. The data points are shown by 
the diamond (SI) and triangle (N2front). The expected colors for both isothermal 
f §5,ip and power-law-thermal ( ^5.2j) cases are indicated as open circles (o) and filled 
circles (•), respectively. OPR=3.0 is assumed for both cases. The different types of 
black lines connect points of equal n(H2) and the LTE case. The grey solid lines 
connect points of equal power-law index (b) or equal temperature (T). The values 
for the power-law index and temperature are also indicated. 
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Fig. 5. The IRC color-color diagrams for the cloud SI with the expected colors for 
isothermal cases of various OPRs (cf. £15.ip . The rest of the graph is the same with 
Figure 0J 
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Fig. 6. Fractional contribution of 12 CO v = 1 — emission lines to the IRC N4 band 
as a function of n(H2). Filled circles and open circles are for b = 4.0 and b = 5.0, 
respectively. Xh is the fractional abundance of atomic hydrogen to molecular hy- 
drogen, iV(H I)/7V(H2). The fractional abundance of CO to H2 is assumed to be 
10~ 4 . Overall, in the range of n(H 2 )=10 3 - 10 6 cm , the fractional contribution of 
12 CO v = 1 - to the IRC N4 band is less than 0.1. 
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Fig. 7. (top-left) The population diagram of the cloud SI, derived from the IRC 
color-color diagram (Fig. H|) using the power-law admixture model (cf. § I5.2D . The 
vibrational levels of v = 0,1,2 are designated by circles, triangles, and squares, 
respectively. The circles show the "ankle-like" curvature, which turn slightly up- 
ward around Ej ~ 3000 K as Ej increases. The filled circles are the levels which 
contribute to the corresponding IRC bands (cf. Table [4]) . The filled triangle is the 
level of (f, J)=(l,3), the upper level of the emission line H2 v = 1 — > S(l). 
The grey filled triangle with an error bar is the population of (v , J)=(l,3) derived 
from the observed H2 v = 1 — » S(l) intensity, extinction corrected, (top-right) 
Two temperature LTE fitting of the H2 levels of the cloud SI, which contribute 
to the corresponding IRC bands (the filled circles of the left figure). OPR=3.0 is 
adopted. The obtained fitting parameters are (T w , N W )=(A61 K, 7.1 x 10 19 cm -2 ) 
and (Th, Nh)=(1225 K, 1.6 x 10 18 cm -2 ). See the section [6.3.1l for a further descrip- 
tion, (bottom-panels) The same plots for the N2front in Cloud N (Paper I). 
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Fig. 8. The expected IRC col o rs fo r a partially dissociative J-shock model 
(iBrand et all ll98Sl : iBurton et all Il989h . The axes represent the ratio of the in- 
tensities in the corresponding IRC bands. The data point is shown by the diamond. 
The connected open-circles have the same fractional CO abundance to N(H.2), X C o] 
the four lines correspond to X co = 0, 10 -6 , 10 -5 , 10~ 4 , from left to right. The pres- 



sure increases from 
line. 



10° cm K to 10 cm K, by a factor of 10 , along each 
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Fig. 9. (top) The schematic description for the bow shock picture (bottom) The 
schematic description for the shocked clumpy ISM picture 
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Fig. 10. The IRC color-color diagram for N2front and Orion Molecular Cloud-1: 
(leftmost) N4/S7 versus S7/S11, (middle) N4/S7 versus S11/L15, and (rightmost) 
S7/S11 versus S11/L15. The grids are the expected colors from the thermal admix- 
ture model (cf. section [5^2]) . Depending on which IRC bands are used for the col- 
or-color diagram, the obtained parameters, b and n(H2), are different for each data 
point; the shorter-wavelength bands return the higher-6 and higher-ra(H2) (cf. N4/S7 
vs. S7/S11 and S7/S11 vs. S11/L15 diagrams). See section 16.3.21 for detail. 
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Table 1 



Summary of the AKARI IRC Observations 


Channel 


Filter 


Wavelength 


Imaging 


Data ID 






coverage" 


Resolution (r) 




(pixel size) 




w 


(FWHM, ") 




NIR 


N3 


2.7-3.8 


4.0 


1402804 


(1.46" x 1.46") 


N4 


3.6-5.3 


4.2 


1402804 


MIR-S 


S7 


5.9-8.4 


5.1 


1402804 


(2.34" x 2.34") 


Sll 


8.5-13.1 


4.8 


1402804 



a Defined as where the res ponsivity is l a rger t han 1/e of the peak 
for the imaging mode. See lOnaka et al.l (|2007l ). 



Table 2 

Results toward Cloud S 



Region 


N4 


S7 


Sll 


N4/S7 


S7/S11 


H 2 v = 1 -> S(l) a 




(MJy sr" 1 ) 


(MJy sr" 1 ) 


(MJy sr" 1 ) 






(erg s — 1 cm -2 sr - 1 ) 


SI 


0.16±0.02 


1.05±0.01 


1.56±0.03 


0.15±0.02 


0.67±0.01 


(5.9±0.2)xl(T 6 



a Extinction-corrected intensity with JV(H)=3.5 X 10 21 cm 2 (Ay = 1.8 mag for Ry=3.X). See text for 
detail. 



Table 3 

Derived Parameters for the Power-law Admixture Model" and the Predicted H2 
v = 1 — ► S(l) Intensity 

Region n(H 2 ) b N(H 2 ; T > 100 K) predicted H 2 v = 1 -> S(l) 

(cm -3 ) (cm -2 ) (erg s _1 cm -2 sr -1 ) 

SI (3.9l 2 ;^)xl0 4 4.2l°;J (2.8^J)xl0 21 (1.5±g;l) x 10" 6 
"See section 15.21 for the detailed model description. 
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Table 4 



Derived Contribution of H2 line emission to the IRC bands 




\ A/n "\ 1 c*y\ o"\~ n 

V V dv ClCllg Lll 


\j jjjjci kj late ±-/iici i^y 




V V LAgllL 


/(J KJKJILvL 1UU 






{X\) 








Ho 77 — n — n 1 ^ 

jri-2 u — ^ yj \± ± j 


A 1 SI 


±0 1 uo 
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^fi9 
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xj 7 , — n n 

11 2 £/ — u — U O^IUJ 


4 41 n 


1 1 QzLD 
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U.O / £i 


R 
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12pf) 7 , _ 1 n 

v v/ V — 1 — U 




ouou 




U.OjU 




Ho a — — 

112 1/ — U VJOIe/l 






XN 4: 


u.000 


49 


H 2 u = - 5(8) 


5.053 


8677 


N4 


0.285 


19 


H 2 v = - 5(7) 


5.511 


7197 


N4 


0.070 


24 


H 2 u = - 5(6) 


6.109 


5830 


S7 


0.346 


8 


H 2 u = - 5(5) 


6.909 


4586 


S7 


0.530 


52 


H 2 v = - 5(4) 


8.026 


3474 


S7 


0.961 


39 


H 2 u = - 5(3) 


9.665 


2504 


Sll 


0.921 


79 


H 2 v = - 5(2) 


12.279 


1682 


Sll 


0.610 


20 



a In units of 10 4 MJy sr 1 /(erg s 1 cm 2 sr 1 ). See the text for the description. 



The contributions are given to the nearest integer. Hence, their sum can be less 
than 100%. 



Table 5 

New Intensity Determination for N2front in Cloud N 



Region 


N4 


S7 


Sll 


L15 




(MJy sr" 1 ) 


(MJy sr" 1 ) 


(MJy sr" 1 ) 


(MJy sr" 1 ) 


N2front 


0.108±0.003 


0.67±0.02 


0.89±0.04 


0.50±0.09 
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